Solving the Cooling Flow Problem of Galaxy Clusters 
by Dark Matter Neutralino Annihilation 
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Recent X-ray observations revealed that strong cooling flow of intracluster gas is not present in 
galaxy clusters, even though predicted theoretically if there is no additional heating source. I show 
that relativistic particles produced by dark matter neutralino annihilation in cluster cores provide a 
sufficient heating source to suppress the cooling flow, under reasonable astrophysical circumstances 
including adiabatic growth of central density profile, with appropriate particle physics parameters 
for dark matter neutralinos. In contrast to other astrophysical heat sources such as AGNs, this 
process is a steady and stable feedback over cosmological time scales after turned on. 

PACS numbers: 95.35.+d, 98.65.Cw 



Diffuse thermal X-ray emission by bremsstrahlung of 
intracluster gas at a temperature of ^ 10 keV has been 
observed from galaxy clusters for many tens of years. 
More than half of clusters are called cooling flow (CF) 
clusters, since cooling time of central cores is less than 
the Hubble time, and theorists predicted the existence of 
strong CF in such systems with a rate of > IQQMq yr~^. 
However, recent X-ray observations failed to reveal evi- 
dence of CFs, requiring that, somewhat ironically, there 
must be some heat source to suppress CFs in "cooling 
flow clusters" Q. Required amount of heating is ^ 10^^ 
erg/s over a time scale of the cluster age ('^ 10^° yr). 

Thermal conduction is probably playing a role, espe- 
cially by preventing the gas from becoming thermally un- 
stable, but a fine tuning is necessary and thermal con- 
duction alone does not successfully explain all clusters^- 
Another heating source popularly discussed is active 
galactic nuclei (AGNs)[3, but efficiency must be very 
high (> 10% of the black hole rest mass energy) 0. Gen- 
erally AGNs are intermittent activity, and accretion rate 
is likely determined by dynamics of small region around 
supermassive black holes (SMBHs). Hence it might be 
somewhat surprising if all clusters are kept stable over > 
100 kpc scale by feedback of central AGNs. 

Clusters are gravitationally dominated by the cold 
dark matter, for which the leading candidate is the light- 
est supersymmetric (SUSY) particles, plausibly the neu- 
tralino X- The neutralino mass is limited in the range 
30 GeV < < 10 TeV, and thermally averaged an- 
nihilation cross section is related to the relic density as 
{<7v) ~ 3 X l{)-'^y{n^h'^) cm-3s-i (see, e.g., for a re- 
view). Detectability of annihilation products from high 
density regions such as the Galactic center (GC) has been 
widely discussed (|^ Q and references therein). Here 
I consider a possibility that annihilation products may 
contribute to heating of intracluster gas. Exotic particle 
dark matter interacting with baryons has been proposed 
to solve the CF problem^, but our scenario is based 
on theoretically better-motivated neutralino dark mat- 
ter. Neutralino annihilation in galaxy clusters has been 



considered by ref.|3l, to explain diffuse radio halos ob- 
served in some clusters. Correlation of radio halos with 
merging clusters 9], however, indicates that such halos 
are formed by cosmic-ray electrons produced by merger 
shocks. The change of central deiisity profile by SMBH 
was not taken into account in ref. 8], without which the 
annihilation luminosity is too small to solve the CF prob- 
lem. 

Since we are interested in relatively central region of 
a cluster, I use the following form of the dark mat- 
ter density profile at r < r^: p = pQ{r/ro)^'^ with 
To — 0.5 Mpc and po = 10^^^ g cm^'^, for a typ- 
ical rich cluster of M15 = Mci/{10^^Mq) - 1. In 
latest numerical simulations of cosmological structure 
formation, the density profile around the center gen- 
erally becomes cusps with 7 ~ 1-1.5^3- Here I use 
7=1, since 7 = 1.5 is not supported by a recent X-ray 
observation |ll|. It is easy to see that neutralino anni- 
hilation has only negligible effects on cluster energetics, 
when simply this density profile is applied. The con- 
tribution from the cusp to annihilation rate (oc p^) is 
convergent with r — > and unimportant when 7 < 1.5. 
Therefore it is a good approximation to estimate the to- 
tal annihilation luminosity, L-^^, using a mean cluster 
density, po- It becomes L^^ ~ 2{av) poMdC^ /m^ 2 x 
10^" {a v) ^261^2^ M 15 erg/s, where TO2 = ^^/(lOOGeV) 
and (cru)_26 = (cru)/(10~^^cm'^s~-^). This is about five 
orders of magnitude lower than that required to suppress 
the CF. 

However, the situation drastically changes if there is 
a density "spike" associated with a central SMBH. Adi- 
abatic growth of a SMBH at the center of a preexist- 
ing halo produces a spike in the density profile within 
r < Ts from the center, which is even steeper than orig- 
inal CUSPS [1^. The power-law index of spike density 
profile, 7s, is related to that of the original cusp as: 
7s = (9-27)7(4-7) = 2.33-2.4 for 7 = 1-1.5. The spike 
radius is a radius within which the enclosed mass of 
the halo is the same as the central SMBH. Since 7s > 1.5, 
the volume integration of the annihilation rate is diver- 
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gent with r ^ 0, and hence a huge enhancement of L-^^. 
is possible from the very central part of the halo. 

If this enhancement is happening in the center of our 
Galaxy, much stronger flux of annihilation gamma-rays 
and cosmic-rays are expected than previously thought, 
and already existing experimental/observational limits 
exclude a considerable SUSY parameter spacejl^- How- 
ever, several authors have argued that such a process 
is rather unlikely to occur in the GC(13i. <14ii |. The ini- 
tial SMBH mass before adiabatic growth should be much 
smaller than the final one, and it must be placed to the 
dynamical center of the halo, but the dynamical friction 
time for that may be larger than the bubble time. The 
mass density of the GC is dominated by baryons rather 
than by dark matter, and violent processes such as star 
formation and supernova explosions might disrupt cold 
orbits of dark matter particles required for the spike for- 
mation. Mergers between halos containing SMBHs may 
lead to a flatter density profile than original cusps. 

However, the center of galaxy clusters, especially those 
having cooling cores within the Hubble time, appear to 
be the best site for the spike growth to happen. All CF 
clusters seem dynamically well evolved systems and a sin- 
gle giant cD galaxy is placed on the grav itational center 
where X-ray surface brightness peaks [l5j- Recent Chan- 
dra observations (e.g., [ll|) have shown that the dark 
matter dominates the mass density with a profile consis- 
tent with 7=1 down to kpc, and probably further 
down to the SMBH scale. If strong CF with the the- 
oretically predicted accretion rate (~ lO^M© yr~-^) oc- 
curred at some epoch in early cluster evolution and this 
accretion is used to feed the SMBH of a cD galaxy, the 
SMBH mass could grow to M.,io = M,/(10^"Mq) - 1 
within > lO^yr depending on the efficiency of mass ac- 
cretion on to the SMBH. This is much shorter than typ- 
ical cluster age, while the orbital period at the spike ra- 
dius rs — LbM^^Q kpc is 5.7 x 10'' Af^^^^ yr, which is 
shorter than the SMBH growth time scale and hence it 
satisfies a requirement for the adiabatic growth. The fi- 
nal SMBH mass is much larger than initial mass that is 
probably similar to typical SMBH mass in normal galax- 
ies {10^~^M(7)), and hence another condition of adiabatic 
growth mentioned by ^3 is satisfied. 

Thus, it seems reasonable to suppose that the adia- 
batic spike growth occurred in the past, within ^ kpc 
of the cluster center where even the latest X-ray satel- 
lites cannot resolve the density profile. Dominant en- 
ergy production by annihilation occurs in the central 
core, where the core density pc is limited by annihila- 
tion itself over typical cluster age, td = 10^°iio yr, as: 
Pc{crv)tci/m^ ~ 1. Equating this core density and the 
spike density profile, p — PQ{r/rs)~'^' (rs/ro)^"', I find the 
core radius Tc — O.ITM^^^q m^'^^'^ {(jv)'^_IIq t^^^ pc. (It is 
much larger than the Schwarzschild radius of the SMBH, 
0.95Af,_io mpc.) The annihilation luminosity within this 



core radius is given by: 

= 1.9 X UfHf.i;„m-''\a^)'!l,t-^'' erg s"'. 

This is a rate within rc, and adding integration at r > Tc 
increases the rate by a factor of 2.8. The above equation 
suggests a modest time evolution of L^^^ oc and 
time average over the cluster age is increased by a factor 
of 7/2. Then I finally obtain L^-^ ~ 10"*^ erg/s, which 
is very close to the number required to heat the cooling 
cluster cores. In contrast to other astrophysical heating 
sources, this process is stable with only mild time evolu- 
tion once it is turned on, unless the spike is disrupted by 
violent events such as major mergers. When such events 
destroy the spike, the cooling core may also be destroyed, 
changing a CF cluster into the other category of non-CF 
ones. 

It must be examined whether this energy production is 
efficiently converted to the heat of intracluster gas. The 
annihilation products are eventually converted to stable 
particles. I used the DarkSUSY package^ to calcu- 
late the amount and energy spectrum of these annihila- 
tion yields, and find that, with only weak dependence on 
SUSY parameters, about 1/4, 1/6 and 1/15 of the to- 
tal annihilation energy goes to continuum gamma-rays, 
e^'s, and pp's, respectively. The other energy goes to 
neutrinos, which are not useful for heating. The spectral 
energy distribution of particles per logarithmic interval, 
e^dN/dt, peaks at about 0.05, 0.05, and 0.1 times m^c^, 
in the same order. Therefore most of the annihilation 
energy will be carried away by particles of about 5 GcV, 
for example, for ^ 100 GeV. 

First I consider the fate of electrons/positrons. For 
simplicity, I assume that a fraction /± ~ 1 of the to- 
tal annihilation energy is given to and all e^'s have 
the same energy of eo = e±/(l GeV) ~ 1. They are 
produced at very dense environment, and their relativis- 
tic motion would result in much higher pressure than 
the environment. Therefore e^'s would expand until 
their pressure becomes comparable with the intraclus- 
ter pressure. Buoyancy may result in intermittent for- 
mation of bubbles of relativistic particles. The den- 
sity in the bubbles can be written by external pressure, 
as: n± ^ -Pcxt/e± 6.3 x 10~^P_9 ep^cm"'^, where 
P_9 = Pcxt/(10^^ erg cm~^) ~ 1. This value and gas 
density n_i = n/(0.1 cm~'^ ) ^ 1 are taken from observed 
values within r < 1-10 kpcjuj. 

Heating of intracluster gas by cosmic ray electrons pro- 
duced by AGNs has been discussed in literature [tR Hif . 
and similar treatments can be applied to estimate the 
energy loss time scale of e^'s, t±. The upper limit on t± 
is given by the ordinary Coulomb collisions with back- 
ground gas: r±.cc = 5.1 x lO^nlJeo yr. This should be 
compared with radiative energy loss time scale by inverse- 
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Compton scattering (ICS) of the cosmic microwave back- 
ground (CMB) photons: T±^ic = 1.2 x IO^Cq^ yr, and syn- 
chrotron radiation: tj- sync = '''±,ic (B/Bcmb) ^, where 
^CMB = 3.30 /zG. I also found that stellar radiation en- 
ergy density in cD galaxies is comparable to that of CMB, 
estimating it by Ust ~ I/cD/(47rr^Qc), where Lcd and Vcu 
are typical observed stellar luminosity and spatial size, 
respectively P^. Therefore ICS of stellar photons should 
also have a comparable effect. Bremsstrahlung loss time 
scale is r±_br ~ 5.7 x lO^nlJ yr and annihilation time 
scale of positrons is ~ 1.0 x lO^^nlJeo yr. Comparing 
these time scales, it can be seen that a considerable part 
of energy can be converted into thermal energy. 

Furthermore, even more efficient energy loss is pos- 
sible when collective effects of plasma, such as the rel- 
ativistic two-stream instability, are important In- 
jection of huge amount of relativistic e^'s within the 
small core radius Tc would lead to strong wind of these 
particles, making their momentum distribution strongly 
anisotropic, which is necessary for the collective effect. 
Using formulae given in these references, I found rj-.tsi — 
3.1 x 103(n±/n)-2(„/cm~3)-o.5g2 g = 8.3 x 10^PZ^nhl4 

yr. Dependence on e± is large, and almost all energy 
could be quickly converted to thermal energy when e± < 
1 GeV. 

The annihilation energy given to pp^s would also be 
partially converted into thermal energy by at least the 
rate of Coulomb collisions that is similar to t±^cc- Some 
fraction of energy may be lost by inelastic interaction 
with ambient thermal protons, with a time scale of Tpp ~ 
(ncTppc)^^ = 3.3 X lO^nlJ yr, and secondary e^'s pro- 
duced by this interaction would be again used for heat- 
ing. It is likely that primary annihilation gamma-rays 
do not contribute much to heat the cooling gas, because 
the optical depth to Compton scattering is ~ 10~^ <C 1. 
However, if there are dense clouds having large optical 
depth around the annihilation core, the Compton heat- 
ing may also have considerable effect, as considered for 
heating by AGNs (e.g., [H). 

What is the preferred value of the neutralino mass in 
this context? The peak of e^dN/de should be less than ^ 
a few GeV, otherwise the energy loss time scale to heating 
(r-|-_cc or r±^tsi) becomes longer than that for radiative 
loss, leading to inefficient heating. Combined with the 
relation L^^ oc m^^^^, the neutralino mass should be < 
100 GeV for the proposed process to efficiently work. 

Now I discuss the observability of any signature of 
the neutralino annihilation. Typically about ~30 con- 
tinuum gamma-rays are produced at > 100 MeV per 
annihilation^^, and expected gamma-ray flux becomes 
F^(> lOOMeV) ~ 7 X 10-^^45 m^^d^^ cm-^s"!, where 
L45 = L^x/iW^^erg s-i) and da = d/(100Mpc) is a 
typical distance to nearby rich clusters. This is just al- 
most the same as the EGRET sensitivity limit, and hence 
the prediction is marginally consistent with no reported 



gamma-ray detection from nearby galaxy clusters |20j. 
In fact, there are many positional coincidences between 
known galaxy clusters and unidentified EGRET sources, 
and the detection from clusters has not yet been claimed 
because of low statistical significance. The next gener- 
ation gamma-ray satellite, GLAST will very likely de- 
tect the continuum gamma-ray flux as steady and point 
sources at cluster centers. 

Even if the continuum gamma-rays may be discrimi- 
nated from other astrophysical sources by spectrum, vari- 
ability, and/or extension, a conclusive evidence would 
come from line gamma-rays. Line flux (in cm^^s^^) is 
expected to be smaller than the continuum flux by a fac- 
tor of about 15 X 10^, where the former factor (30/2) 
comes from the ratio of photon number produced per 
annihilation and the latter comes from the branching ra- 
tio into XX ~^ 77 or Zj modes 0. Following the line 
sensitivity estimate given in this reference, I found that 
the line flux expected at < 100 GeV may produce sev- 
eral photons from a cluster center for flve-year operation 
of GLAST, compared with negligibly small background 
rate of ^ 10~^ events within angular and energy resolu- 
tion. Co-added analysis of many CF clusters would even 
increase the sensitivity. Future air Cerenkov telescopes 
may have even better sensitivity for the line flux, but the 
threshold energy must be lower than ~ 100 GeV since we 
expect < 100 GeV. 

Since the flux enhancement by the density spike is so 
drastic, the best target to search neutralino annihilation 
may be nearby CF clusters, rather than the GC or nearby 
galaxies. On the other hand, L^^ without the density 
spike (^ 10^° erg/s) is below the sensitivity limit even 
for the next generation gamma-ray telescopes. 

Equally intriguing is the detectability of synchrotron 
radiation in the radio bands, by pairs produced by 
annihilation. It is known that central cD galaxies in 
CF clusters have a higher probability of becoming radio 
sources than those in non-CF clusters, and there seems 
a correlation between radio power and CF strength [ij , 
suggesting that the annihilation may be responsible 
for radio emission. Typical radio luminosity from CF 
cluster cores is ~ 

2041-42 erg/sHH, which is about 10^ 
times lower than the cluster X-ray luminosity or L-^^. 
This factor can be explained by a few effects as follows. 
The magnetic field strength is poorly known, which may 
be weakened by relativistic bubble formation. Hard X- 
ray tails observed in some cluster spectra are popularly 
interpreted as ICS of CMB photons by cosmic-ray elec- 
trons, and their luminosity is a few hundreds times larger 
than diffuse radio halo luminosity, indicating that syn- 
chrotron is inefficient process compared with ICS 2^. 
The energy loss by heating can be much more efficient 
than radiative loss by the collective effect, and then the 
total radio flux will be further reduced. Characteristic 
synchrotron frequency is Vgync ~ .01le^{B / jiG) GHz. 
Since the injection of relativistic occurs into relatively 
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narrow energy range around ~ Q.OBm^c^ compared with 
broad power-law spectrum of cosmic-ray electrons, most 
of radio emission may occur at < GHz, which is out of 
typical observing frequencies. 

Here I give a size estimate of the radio emitting region. 
It is expected that the size of bubble, r?,, is determined 
by the pressure balance and the lifetime of these parti- 
cles, as: rt = (3T±iV±/47rn±)-^/'^, where N± — f±L^^/e± 
is the production rate. Using r±_tsi for the lifetime 
of e±'s, I found n = 13 n]/^eo^^Pl^ fl/^Ll^^ kpc. Typ- 
ical extension of the radio emission associated with cD 
galaxies is ^ 5-10 kpc 21]. Large scale radio halos (> 
Mpc) are found only in non-CF clusters like Coma, but 
CF clusters often have less extended "minihalos" (< 100 
kpc) in which a strong radio cD galaxy is centered [23| . 
It may be speculated that the bubble corresponds 
to the central strong radio source, while radio miniha- 
los are made by leaking from the bubble, perhaps by 
diffusion. The morphology of such radio galaxies and 
minihalos in CF clusters is poorly coUimated and more 
spherical, compared with clear bi pol ar jet-like structures 
generally found in radio galaxies |2l|. Such a trend may 
be difficult to explain if CF is suppressed by AGN jets, 
while it is naturally understood if the steady, isotropic 
energy production by neutralino annihilation is responsi- 
ble. If this interpretation is true, it means that we have 
already observed, though indirectly, the dark matter for 
many tens of years! 

Annihilation pairs would also produce a similar lu- 
minosity to synchrotron radiation in the X-ray and MeV 
bands by ICS of CMB and stellar photons, respectively, 
but they are difficult to detect because of the strong ther- 
mal X-ray emission and difficulty of MeV observations. 

The simple estimate of L^^ without the density spike, 
~ 10""^ erg/s, is comparable to the power of diffuse radio 
halos in merging clusters and hence neutralinos may be 
the origin, as pointed out by ref.[l|- A necessary condi- 
tion for this is that the synchrotron must be the dominant 
energy loss process. However, as argued above (and refs. 
H^), observed hard X-ray tails in cluster spectra indicate 
that synchrotron is rather inefficient process, and diffuse 
synchrotron halo by annihilation e^'s is likely still under 
the current detection limit. 

It is expected that stars in cD galaxies are also affected 
by the adiabatic growth, and it might be interesting to 
seek for any signature in the central luminosity density 
profile of cD galaxies by high resolution optical obser- 
vations. In fact, adiabatic growth by SMBHs is one of 
the proposed explanations of cusps seen in central sur- 
face brightness profile of elliptical galaxiesjl^. However, 
there is no reason to believe that stellar density profile 



should be the same as that of dark matter, and large ellip- 
tical galaxies generally have flat stellar density cores |24j. 
A stellar density spike should be weak and may not be 
detectable, if it is formed from such a flat density core. 
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